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polynomials in TT ^ and into the formal enveloping series 
with rational coefficients only 
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Abstract 

In this work, two new series expansions for generalized Euler's constants (Stieltjes constants) 7 m are 
obtained. The first expansion involves Stirling numbers of the first kind, contains polynomials in ti~^ 
with rational coefficients and converges slightly better than Euler's series y The second expan¬ 
sion is a semi-convergent series with rational coefficients only. This expansion is particularly simple 
and involves Bernoulli numbers with a non-linear combination of generalized harmonic numbers. 
It also permits to derive an interesting estimation for generalized Euler's constants, which is more 
accurate than several well-known estimations. Finally, in Appendix A, the reader will also find two 
simple integral definitions for the Stirling numbers of the first kind, as well an upper bound for them. 
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I. Introduction and notations 

1.1. Introduction 

The ^-function, which is usually introduced via one of the following series. 
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is of fundamental and long-standing importance in modern analysis, number theory, theory of special 
functions and in a variety other fields. It is well known that ^(s) is meromorphic on the entire complex 
s-plane and that it has one simple pole at s = 1 with residue 1. Its expansion in the Laurent series in 
a neighbourhood of s = 1 is usually written the following form 
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where coefficients 7 m, appearing in the regular part of expansion ||2]l, are called generalized Euler's 


constants or Stieltjes constants, both names being in useQj Series l|2]l is the standard definition for 7 m 


Alternatively, these constants may be also defined via fhe following limif 
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The equivalence between definitions (|2ll and (|3]l was demonstrated by various authors, including 
Adolf Pilz Thomas Stieltjes, Charles Hermite ill vol. I, letter 71 and following], Johan Jensen 


|87tl89l], Jerome Franel llSbIl . 


orgen P. Gram Godfrey H. Hardy I. Srinivasa Ramanujan j^j. 


William E. Briggs, S. Chowla J 24 j and many others, see e.g. lllbl.ll7bLl84l,ll28ll . It is well known that 70 = 


7 Euler's constant, see e.g. 111281] , il9l. Eq. (14)]. Higher generalized Euler's constants are not known 
to be related to the "standard" mathematical constants, nor to the "classic" fimctions of analysis. 

In our recent work lll 8 ll . we obtained two interesting series representations for the logarithm of the 
T-function containing Stirling numbers of the first kind Si (n, k) 
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as well as their analogs for the polygamma fimctions The present paper is a continuation of 

this previous work, in which we show that the use of a similar technique permits to derive two new 


^The definition of Stieltjes constants accordingly to formula (2) is due to Godfrey H. Hardy. Definitions, introduced by 
Thomas Stieltjes and Charles Hermite between 1882-1884, did not contain coefficients ( — 1)'" and m! In fact, use of these 
factors is not well justified; notwithstanding, Hardy's form (2) is largel y ac cepted and is more frequently encountered in 
modern literature. For more details, see a vol. I, letter 71 and following], 11 id p. 562], il^, pp. 538-539]. 

®Some authors use the name generalized Eu ler’s constants for other constants, which were conceptually introduced and stud¬ 
ied by Briggs in 1961 and Leh merin 19 75 jll4i]. They were subsequently rediscovered in various (usually slightly different) 

forms by several authors, see e.g. ll73l,ll4dll9Cll . Further generalization of both, generalized Euler's constants defined accord¬ 
ingly to (2} and generalized Euler's constants introduced by Briggs and Lehmer, was done by Dilcher in li^ . 

^Both series converge in a part of the right half-plane S Fig. 2] at the same rate as E [n ln’"n) ^, where m = 1 for InT(z) 
and Y(z), m = 2 for Yi(z) and Y 2 (z), m = 3 for Y 3 (z) and Y 4 (z), etc. 
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series expansions for generalized Euler's constants 'Ym, both series involving Stirling numbers of the 
first kind. The first series is convergent and contains polynomials in tz~^ with rational coefficients (the 
latter involves Stirling numbers of the first kind). From this series, by a formal procedure, we deduce 
the second expansion, which is semi-convergent and contains rational terms only. This expansion is 
particularly simple and involves only Bernoulli numbers and a non-linear combination of general¬ 
ized harmonic numbers. Convergence analysis of discovered series shows that the former converges 
slightly better than Euler's series in a rough approximation at the same rate as 


“ In'”Inn 
^3 n^ln^n ' 


m = 0,1,2,... 


The latter series diverges very quickly, approximately as 
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1.2. Notations and some definitions 

Throughout the manuscript, following abbreviated notations are used: 7 = 0.5772156649... for 
Euler's constant, 7 m for mth generalized Euler's constant (Stieltjes constant) accordingly to their def¬ 
inition ||2]lj^ (^) denotes the binomial coefficient C[h B„ stands for the nth Bernoulli numberj^ H„ and 

(s) 

H„ ' denote the nth harmonic number and the nth generalized harmonic number of order s 
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respectively. Writings [xj stands for the integer part of x, tgz for the tangent of z, ctgz for the 
cotangent of z, ch z for the hyperbolic cosine of z, sh z for the hyperbolic sine of z, thz for the hy¬ 
perbolic tangent of z, cthz for the hyperbolic cotangent of z. In order to avoid any confusion be¬ 
tween compositional inverse and multiplicative inverse, inverse trigonometric and hyperbolic func¬ 
tions are denoted as arccos, arcsin, arctg,... and not as cos“^, sin“^, tg“^,_ Writings r(z) and 

^(z) denote respectively the gamma and the zeta functions of argument z. The Pochhammer sym¬ 
bol (z)n, which is also known as the generalized factorial fimction, is defined as the rising factorial 
(z)„ = z(z-h l)(z-k 2) ■ ■ ■ (z-k n - 1) = r(z-k n)/r(z)13§ For sufficiently large n, not necessarily 
integer, the latter can be given by this useful approximation 

, , r 6z^ —6z-kl Sbz"^ — 120z^-k 120z^ — 36z-k 1 

= r(z)c>^ 12n +-288^^^-+ ) 

( 6 ) 

n^-T{n) f z(z-l) z(z - l)(z - 2)(3z - 1) 

“ r(z) i ^ 2 n + 24n2 ^ ^ 


5ln particular 71 = -0.07281584548..., 72 = -0.009690363192..., 73 = +0.002053834420 .... 

^In particular Bq = +1, Bi = - j, B 2 = + 5 , B 3 = 0, B 4 = B 5 = 0, Be = +^, B/ = 0, Bg = -^, B 9 = 0, Bio = +^/ 
Bii = 0, Bi 2 = see Tab. 23.2, p. 810], jl09l. p. 5] or p. 258] for further values. Note also that some authors 

may use slightly different definitions for the Bernoulli numbers, see e.g. 1^, p. 91], Ei pp. 32, 71], izil p. 19, n° 138] or El 
pp. 3-6]. 

^For nonpositive and complex n, only the latter definition (z)„ = r(z + n) /r(z) holds. 

^Not e tha t some writers (mostly G erma n-speaking) call such a function/flcu/fe analytiqiie or Facultat, see e.g. Il 57 ll , EE 
p. 1 86], I159I. vol. II, p. 12 ], 1^, p. 119], liodl . Other names and notations for (z)„ are briefly discussed in 1^, pp. 45^7] and in 
S pp. 47-48]. 
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which follows from the Stirling formula for the F-functionl^ Unsigned (or signless) and signed Stirling 
numbers of the first kind, which are also known as factorial coefficients, are denoted as |Si (n,l) \ and 
Si{n,l) respectively (the latter are related to the former as Si{n,l) = ( — l)”^^|Si(n,Z)|)113 Because 
in literature various names, notations and definitions were adopted for the Stirling numbers of the 
first kind, we specify that we use exactly the same definitions and notation as in Section 2.1], 
that is to say \Si{n,l)\ and Si{n,l) are defined as the coefficients in the expansion of rising/falling 
factorial 


n-l n 00 
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k^O ^ /=i /=o 
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(7a,b) 


respectively, where z G C and n ^ 1. Note that if I [1, n], where I is supposed to be normegative, 
then Si (m, Z) = 0, except for Si (0,0) which is set to 1 by convention. Alternatively, the same numbers 
may be equally defined as the coefficients in the following MacLaurin series 
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Signed Stirling numbers of the first kind, as we defined them above, may be also given via the fol¬ 
lowing explicit formula 


„ ( 2 n-Z)!”^' 

= V-utE 


1 k ^_.^y^n-l+k 
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Z G [1, n], which may be useful for the computation of ^(n, Z) when n is not very large A ll three 
above definitions agr ee with those adopted by Jordan Chapt. IV], |90, Riordan il47i p. 70 et 
seq.], Mitrinovic I125ll . Abramowitz & Stegim yl n° 24.1.3, p. 824] and many others (moreover, mod¬ 
ern CAS, such as Maple or Mathematica, also share these definitions; in particular Stirlingl(n,l) 
in the former and StirlingSl [n,l] in the latter correspond to our Si(n, i))0 Kronecker symbol (or 
Kronecker delta) of arguments Z and k is denoted by {Sij,- = 1 if Z = Zc and li; j- = 0 if I k). 
Re z and Im z denote respectively real and imaginary parts of z. Letter i is never used as index and is 
V—1 ■ The writing resz=fl /(z) stands for the residue of the function /(z) at the point z = a. Finally, by 
the relative error between the quantity A and its approximated value B, we mean [A — B) /A. Other 
notations are standard. 


® A simpler variant of the above formula may be found in ilvA . 

^‘^There exist more than 50 notations for the Stirling numbers, see e.g. 1^ , 1^ pp. vri-viri, 142,168], llOlI pp. 410—422], 
Sect. 6.1], and we do not insist on our particular notation, which may seem for certain not properly chosen. 

^^From the above definitions, it follows that: Si(l,l) = +1, Si(2,l) = —1, Si(2,2) = +1, Si(3,l) = +2, Si ( 3,2) = —3, 
Si(3,3) = +1, ... , Si(8,5) = “I960, ... , Si(9,3) = +118124, etc. Note that there is an error in Stirling's treatise 11721] : in the 
last line in the table on p. 11 jl^ the value of | Si (9,3) | = 118124 and not 105 056. This error has been note d by Jacques Binet 
Ei p. 231], Charles Tweedie Il79l, p. 10] and some others (it was also corrected in some translations of jl72ll 'l. 

quick analysis of several alternative names, notations and definitions m ay be found in works of Charles Jordan 1^, 
pp. vii-viii, 1 and Chapt. IV], Gould 1^1^ . and Donald E. Knuth 1^, Sect. 6.1], llOll, pp. 410^22], 
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II. A convergent series representation for generalized Euler's constants 7,,, involving Stirling num¬ 
bers and polynomials in 


III. Derivation of the series expansion 

In 1893 Johan Jensen 1881,18911 by contour integration methods obtained an integral formula for the 
^-function 
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1 1 f (1 — ix) ® —(1 + fx)' 


s -1 2 z 


einx _ I 


■ dx, s 7 ^ 1 


which extends iffll to the entire complex plane except s = 1. Expanding the above formula into the 
Laurent series about s = 1 and performing the term-by-term comparison of the derived expansion 
with the Laurent series (|2| yields the following representation for the mth Stieltjes constant 


00 
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which is due to the Jensen and Franellil Making a change of variable in the latter formula 
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where m = 0 , 1 , 2 , 

Now, in what follows, we will use a number of basic properties of Stirling num bers , w hich ca n be 


found in an amoim t sufficient for t he p resent purpose in th e following literature: 


| 1 ^ 


| 1 M, 


Book I, part IL |52l,ll55lll56 . 
lii p. 129], STchaM. IV], 
422], H, Chant VL E" 


91 


171] 


61 


189 


1 1581. P D. 186-187], lisa vol. II PP - 23-31], |lldl32, 


33 


I 


2 lll^, 


131, H PP- 67-78L S Sect. 6. 11. jlOll. pp. 410- 

_, 1141], Ch^t. 4, §3, n° 196-n°210] ||5~p. 60 et seq.], ilOol]. j un p. 70 et seq.], 
Chapt 8 ], S n° 24.1.3, p. 824], H Sect. 21.5-1, p. 824]^, vol. Ill, p. 257], 


11691 vol . 1], 1511, 13^ ' _ . 

J^pp. 91-m im pp. 2862-28651 ini Chapt. 2 1. y.y.lmEal^. E p . 642], 


mSH m[ id EE 0011 mQliii i6ili64l. 154 


148 


149, 76 


lol 


M- Note 


that many writers discovered these numbers independently, without realizing that they deal with the 


the explicit form, this integral formula was given by Franel in 1895 Is^ (in the above, we corrected the original Franel's 
formula which was not valid for m = 0). However, it was remarked by Jensen that it can be elementary derived from 1101 
obtained two years earlier and it is hard to disagree with him. By the way, it is curious that in works of modern authors, see 
e.g. formula 111! is often attributed to Ainsworth and Howell, who discovered it independently much later @]. 
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Stirling numbers. For this reason, in many sources, these numbers may appear under different names, 
different notations and even slightly different definitions!^ 

Consider the generating equation for the unsigned Stirling numbers of the first kind, formula ||8K). 
This power series is uniformly and absolutely convergent inside the disk |z| < 1. Putting I + m — 1 
instead of I, multiplying both sides by {l)m and summing over I = [1, oo), we obtain for the left side 


^ ^ [-ln(l-z)] 

k (l + m-iy. k (^- 1 )! 
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while the right side of lHh), in virtue of the absolute convergence, becomes 
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Writing in the latter — z for z, and then subtracting one from another yields the following series 


In'^(l-z) _ ln'"(l+z) 
1 — z 1 + z 


2 (-l)'"m! 


“ |Si (2A: + 2, m + 1) I 2k+i 

h " 


(14) 


m = 0,1,2 ,..., which is absolutely and uniformly convergent in the unit disk |z| < 1, and whose 
coefficients grow logarithmically with k 


|Si(2k + 2,m + l)| hTk 


(2fc + l)! 


m\ 


fc —>■ 00 , m = 0 , 1 , 2 ,.. 


(15) 


in virtu e of k nown asymptotics for the Stirling numbers, see e .g. |9ll, p. 261],jl92 


p. 824], 11881 p. 348, Eq. (8)]. Using formulae from p. 217], 11631, p. 1395], 1104 p. 425, Eq. (43)], the 


p. 161], Sn^ 24.1.3, 


Actually, only in the beginning of the XXth century, the name "Stirl ing numbers" appeared in mathematical literature 
(mamly, thanks to Thorvald N. Thiele and Niels Nielsen EiEESi, linil. p. 416]). Other names for these numbers include: 
factorial coefficients, faculty's coefficients {Facultatencoejficienten, coefficients de lafaculte analytique), differences of zero and even differ- 

is) 

ential coefficients of nothing. Moreover, the Stirling numbers are also closely connected to the generalized Bernoulli numbers B„ , 
also known as Bernoulli numbers of higher order, see e.g. p. 129], p. 449], p. 116]; many of their properties may be, 

fs) 

therefore, deduced from those of BU- 
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law for the formation of first coefficients may be also written in a more simple form 


Si(2fc + 2,m + l) 
{2k + l)\ 


1 , 


H2k+1 / 



m = 0 
m = 1 

(16) 

m = 2 
m = 3 


For higher m, values of this coefficient may be similarly reduced to a non-linear combination of the 
generalized harmonic numbers. Since expansion ifTH holds only inside the unit circle, it cannot be 
directly used for the insertion into Jensen-Franel's integral formula dTTb . Flowever, if we put in ifT^ 
z = 2 ^ in(l — u), we obtain for the right part 


2 (-l)'”m! 


00 


■E 
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{2k + l)\ 
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see (8^) 
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= 2f(-l)™m!.E^-E 
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Therefore, for m = 0,1,2,..., we have 


“ (-l)'^|Si(2fc + 2,m + l)| ■ |Si(n,2fc + l)| 

n=l ”■ k=0 


{2ny^+^ 
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' In'” 

\ In(l-w)' 
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In'^' 
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k=0 


(27r)2^+i 


(18) 


which uniformly holds in | m| < 1 and also is valid for = 1 Substituting ifTSll into ifl^ and per¬ 
forming the term-by-term integration from u = 0 to u = 1 yields the following series representation 
for mth generalized Euler's constant 


7m - 2 ^m,o + - E TTTm' E 


1 (-l)^|Si(2fc-h2,m-hl)| ■ |Si(n,2A:-hl)| 


(19) 


n=l " ■ “o (27r)2'^+l 

where m = 0,1,2,.... In particular, for Euler's constant and first Stieltjes constant, we have following 


^®The unit radius of convergence of this series is conditioned by the singularity the most closest to the origin. Such singularity 
is a branch point located at u = 1. Note also that since the series is convergent for i; = 1 as well, in virtue of Abel's theorem on 
power series, it is uniformly convergent everywhere on the disc | w | < 1 — £, where positive parameter e can be made as small 
as we please. 
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I Series for yn-Series for yi -Series for yi 



Figure 1 : Absolute values of relative errors of the series expansion for 70, 71 and 72 given by (T^-j 20 ), loga¬ 
rithmic scale. 


series expansions 


1 ^ “ 1 (-l)^-(2fc + l)!- |Si(n,2fc + l)| 

2 ^ ^ w ■ «! (27 t)2*^+i 

1 J_ J_ __ 

2 2tz^ 8tz^ 18 \7r2 4/1^/ 96 \7 t2 271“^/ 

J_/12_l^ 1 / 60 675 225 \ 

600 \ Tz^ Atz^ 47t® / ”*” 4 320 \ 47t^ ) 


7i 
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n ■ n! 


L2”J 
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k=0 


i-lf-{2k+l)\-H2k+i-\Si{n,2k + l) 
(27r)2^+i 


1 

1 

1 / 

^ 1 11 \ 

3 / 

M 11 \ 

27r2 

8n^ 

18 1 

loo 

1 

IfN 

1 

96 \ 

JJ.2 J 


J_/12_^ 137 \_1 /60 2 475 2055 \ 

600 \7r2 Stt^ 167t^/ 4320\7r2 8n^ 167t^/ 


( 20 ) 


respectively. As one can easily notice, each coefficient of fhese expansions confains polynomials in 
Ti~^ with rational coefficients. The rate of convergence of fhis series, depicfed in Fig. [TJ is relafively 
slow and depends, at least for the moderate number of terms, on m: the greater the order m, the 
slower the convergence. A more accurate description of this dependence, as well as the exact value of 
fhe rate of convergence, both require a detailed convergence analysis of ifT^ . which is performed in 
fhe next section. 


11.2. Convergence analysis of the derived series 

The convergence analysis of series ifT^ consisfs in the study of ifs general term, which is given 
by the finite truncated sum over index k. This sum has only odd terms, and hence, by elementary 
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transformations, may be reduced to that containing both odd and even terms 


E(-i) 

k=0 


^ |Si(2fc + 2,m + l)| ■ |Si(n,2fc + l)| 


lin\ 


= ( — 1 ) 2 ( 2 ^+ 1)“2 


(27t)2*^+i 

^( 2 k+i)—^ I Si ( 2 /c + 1 + 1 , + 1 ) I ■ |Si(M, 2 fc + l)| 
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( 21 ) 


k=0 

_ r_iVi. . |Si(^+ !,?« +1)1 ■ |Si(n,/)| 


= oE[i-(-ir]'(-i) 


■1=1 


{Iny 


where, in the last sum, we changed the summation index by putting I = 2k+ 1. Now, from the second 
integral formula for the unsigned Stirling numbers of the first kind, see llA.4b , it follows that 


■ dz 


|Si(Z + l,m + l)| 

i-ir (Z + l)! / 

i 

l 

1 

{2ny 

2ti (m + 1 )! J 

27TZ 

z^ 
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(-i)h(-i)ia- 


1 ) 1 Si (Z “F 1, + 1 ) 1 

(- 1 )- (z + 1 )! r 

i 

1 in’«+l(i 

(2/1)' 

2n (m -F 1)! / 

2nz_ 

z?- 


■ dz 


|z|=r 


where 0 < r < 1 . Therefore, since (Z + 1)! = ^dx taken from 0 to oo, the last sum in ll 2 Tll 

reduces to the following integral representation 


(_l)m 


Inz J \ 2nz J 


ln™+i(l-z) 


dz 


ElSiKOl 

1=1 


ATz{m +1)! 

(-l)m 

47r(m + l)! j z^ 
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/ ix \ / ix \ 
\27rz/ \ 27rzy 


ln'”+i(l-z) 


dz 


xe ^ dx 


( 22 ) 


ln”'+i(l-z) 




\2nzJ„ V 27tz;J 


ix \ 


xe ^ dx\ dz 


The integral in curly brackets is difficult to evaluate in a closed-form, but at large n, its asymptotical 
value may be readily obtained. 

Function 1/T(z) is analytic on the entire complex z-plane, and hence, can be expanded into the 
MacLaurin series 


T(z) 


2 /7^ 3 ^ k II 

= z-h7Z + (Y - Y 2 - E^ “ft/ |z|<<». 


( 23 ) 
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where 



■ 1 
. 1 ^ 
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(- 1 )^ 
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see e.g. 0, p. 256, n° 6.1.34], 1 188 , pp. 344 & 349], 177]. Using Stirling's approximation for the Pochham- 
mer symbol lUJ, we have for sufficiently large n 


/ ix \ / ix \ 


r(n) n ^ ■ r(n) 


= {n-l)\ 


2tizJ„ 
'ixlnn\ ^ 


r(*) 




2nz J 


(24) 


exp 


Inz 




k=l 


ix 

2nz ) 


exp 


zxlnn\ ^ , V f tx \ 
„ , , , air I -— I 

2nz 




Substituting this approximation into the integral in curly brackets from ll22ll . performing the term- 
by-term integratioio and taking into account that z“®r(s) = f x^~^e~^^dx taken over x e [0, oo), 
yields 


oo 

/ 


/ ix \ / ix \ 

A277Z/n V 27Tzj„_ 


xe ^ dx ^ 




1 + 


ilnn 

2nz 


-k-2 


(n- 1)! 


32 in^z^ — 3 In^ n) 


-(-ir 


n —>■ 00 , 


1 - 


f Inn 
2nz 


-k-2' 


(25) 


(47r2z2 -|- Im n) 

where, at the final stage, we retained only the first significant term corresponding to factor k = ll^ 
Now, if |z| ^ 1 — e~^ ~ 0.63, then the principal branch of | ln'”'''^(l — z) | ^1 independently of m and 
argz. Analogously, one can always find such sufficiently large ng, that for any however small e > 0, 

32 in^z^ (47r^z^ — 3 In^ n) 


< e, no, 


(26) 


(47t2z^ -|- Im n) 

on the circle |z| = 1 — e“^ (for example, if e = 1, then ng = 1222; if e = 0.1, then ng = 38597; 
if e = 0.01, then ng = 33220487; etc. )@ Combining all these results and taking into account that 
|dz| = |z| d argz, we conclude that 


n ■ n! 


LAJ (-1)^1 Si(2 k + 2,m +1)1 ■ |Si(n,2k + l)| 
k=0 


(27r)^*^+i 


< 


(27) 


< 


2(1 -e-i)(m-M)! 


< 


C 


n ^ ng, 


Numerical simulations, see Fig. 12 show that this simple inequality, valid for all m, may provide more 
or less accurate approximation for the general term of (l22ll , and this greatly depends on m. Moreover, 


^^Series (23)-(24) being uniformly convergent. 

^^The second term of this sum, corresponding to /c = 2, is 


396 (4/1^2^ — Ln^ «) in n 

{n-iy. -— = (n-iy.-0 


In^n 


(472^2^ + Im n) 
and hence, may be neglected at large n. 

'®Note that for fixed n, the left-hand side of (26) reaches its maximum when z is imaginary pure. 
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I m = Q-m=l - m = 2 -m = 3 



Figure 2: Relative error between the upper bound and the left-hand side in i27l as a function of n for four 
different orders m, logarithmic scale (the curve in long dashes correspond to m = 1, that in short ones to m = 3). 
Results displayed above correspond to C = 1/ (2/t). 


the joint analysis of Figs. [T] and |2] also indicates that first partial sums of series dl^ may behave quite 
irregularly One of the reasons of such a behaviour is that for 1 ^ n ^ 53, absolute value ll26ll increases, 
and it starts to decrease only after n = 540 Notwithstanding, inequality ll27ll guarantees that in all 
cases, the discovered series for jm given by llT9ll converges for large n not worse than Euler's series 
Ln-^. 

Asymptotics ll25l l. as well as the rates of convergence previously obtained in 11811 . both sug gest 
that the exact rate of convergence of series (IT^ may also involve logarithms. For instance, from 1118 . 
Sect. 3], it straightforwardly follows that the rate of convergence of this series at m = 0 is equal to 
J^{n\nn)~^, see also footnote ID Indeed, if we replace the integral in curly brackets from ll22t by its 
first-order approximation (l25ll . and then, evaluate the sum of corresponding residues at Zi 2 = i 


^ z(47r^2^ — 31n^n) ln(l — z) 

\ j'05_j_ i _ 

1=1 (4/1^2^ -|- In^ n) ^ 


In^ n — 471^ 
8 / 1 ^ (471^ -b In^ 


Stt^ In^ i 


^ z(47r^z^ — 31n^n) ln^(l — z) In^ n ■ In (47 t^- b In^ n)-b ... 2-lnlnn 

) res ^- - —5 - = - ^^~,— 


l=\ ^ (47 z^z^ -b In n)" 


87r2(47T2 -b In^ 


^ z(47r^z^ — 31n^n) ln^(l 

) res-^ 

1=1 (47r2z2 -b In^ n) 


z) 3 In^ n ■ [ In^(27T -b i In n) -b In^(27 t — f In n)] -b ... 
167 t2(47t 2 -b In^ n)^ 

3-ln^lnn 

Sn^ln^n 


the circle |z| = 1 — e absolute value has one of its third-order poles at n = 53.08. Other poles are 

located either below n = 1, e.g. n = ^- 27 /( 1 -/ ) 0.02, or are complex. More precisely, all poles of this expression occur at 

n = [cos (27r(l — cos f) isin (27r(l — cos )smip^ where <p = argz. 
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and so on, we find that 

1 JSi(2fc + 2,m + l)| ■ |Si(w,2/c+l)| ^ 

n ■ nl ^ (27t)2^+^ 

(28) 

8i7r2(-l)™ ^ z(47rV-31n2n)ln'"+\l-z) ^ ^ ,^m+l27^ ln'"lnn 

^ nHni +1)1 f (47r2z2 + ln2n)' ^ ^ 

Iz|=r ^ ^ 

in virtue of the Cauchy residue theorem. Of course, this formula is only a rough approximation, 
because poles Zi 2 belongs to the disc |z| = r < 1 only if 1 ^ n ^ 535, while formula ll25ll is a double 
first-order approximation and is accurate only for large n. Furthermore, residues were also evaluated 
only in the first approximation. Fiowever, obtained expression gives an idea of what the true rate of 
convergence of series (IT^ could be, and it also explains quite well why series for higher generalized 
Euler's constants converge more slowly than those for lower generalized Euler's constants. Moreover, 
this approximation agrees with the fact that ifl^ converges not worse than Euler's series, and also is 
consistent with the previously derived rate of convergence for 7 from Jl^ . which was obtained by 
another method. 


III. Expansion of generalized Euler's constants 7m into the formal series with rational coefficients 


III.l. Introduction 

Expansions into the series with rational coefficients is an interesting and challenging subject. There 
exist many such representations for Euler's constant 7 and first of them date back to the XVIIIth 
century. For instance, from the Stirling series for the digamma function at z = 1, it straightforwardly 
follows that 


^2^+^ = 1 + 1 

2 ' 2k 2{N + 1) 2 12 


1 1 1 1 691 

m E2 “ ^ “ 32760 ■ 


(29) 


where 0 < 0 < 1 and N < oo|^ A more general representation of the same kind may be obtained by 
Euler-MacLaurin summation 


n -I T 

7 = E ^ - ln« - — 

I ^ ir 2n 


k=l 


N 

E 

k=l 


^Ik 


0- 


B2N+2 


2k-n^^ 2(N-hl)-n2N+2' 


(30) 


where O<0<1, N<oo and n is positive integer, see e.g. IZll n°377]. Two above series are semi- 
convergent (or divergent enveloping), i.e. they diverge as N —> 00. The first known convergent series 
representation for Euler's constant with only rational terms, as far as we know, dates back to 1790 
and is due to Gregorio Fontana and Lorenzo Mascheroni 

\G„\ 1 1 1 19 3 863 

2^ 362880 


7 


E- 

n—1 


_ 1 1 1 

~ 2 ^^ 


(31) 


^'^This result should be attributed to both Stirling and De Moivre, who originated Stirling series, s ee Il72l. p. 1 351 a nd 13 
resp ectively (for more information on Stirling ser ies, s e e also lil. part II, Chapter VI, p. 466], 13. p. 33], 1251, p. 3291, Il35l, p . Ill ], 
liii §12-33], S §15-05], lii p. 530], S p. 1], [13, EH §4.1, pp. 293-294], [si pp. 286-288], E n‘>6.1.40-n“ 6.1.41], ii3l l. 
Curiously, Srivastava and Choi I 166 I. p. 6], did not notice the trivial connection between this series and the Stirling series for 
the digamma function and erroneously credited this result to Konrad Knopp, in whose book I3l it appears, with a slightly 
different remainder, on p. 527 (Knopp himself never claimed the authorship of this formula). 
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where rational coefficients Gn, known as Gregory's coefficients^^ are given either via their generating 
function 


or explicitly 


i— 7TGI —t ~ ^ Gnz”, 

ln(l+z) 


|z| < 1, 


(32) 


^ 1 ^ Si(n,/) If, , , 

Gn = — > = — {x-n + l)ndx= -■ 

n! g / + 1 n!i 


B 


(n-l) 


_ Gl,n 

(n —l)n! n\ 


(33) 


This series was first studied by Fontana, who, however, failed to find a constant to which it con verges. 
Mascheroni identified this Fontana's constant and showed that it equals Euler's constant jligi pp. 21- 
231- This series was subsequently rediscovered man y tim es, in particular, by Ernst Schroder in 1879 
1 161 , p. 115, Eq. (25a)], W Niels E. Norlund in 1923 |135, p. 244], by Jan C. Klu3wer in 1924 j^ , by 
Charles Jordan in 1929 p. 148], by Kenter in 1999 193], by Victor Kowalenko in 2008 il04 , 103 j. An 
expansion of a similar nature 


7 = 1- E 


C2,n 


251 


19 


19087 


_ 1 5 1 

n^i 14400 1728 2540160 


(34) 


where rational numbers C 2 ,n, known as Cauchy numbers of the second kirn 

z ^ ~ (-lyv^ 

(l+z)ln(l + z) n! 


C2.n=t^^= fix)ndx= |bW 

1=1 ' + ^ i 


(35) 


follows from a little-known series for the digamma function given by Jacq ues B inet in 1839 [I 


p. 257, Eq. (81)] and rediscovered later by Niels E. Norlund in his monograph 11351, p. 244]|f3 Series 


These coefficients are also called (reciprocal) logarithmic numbers, Bernoulli numbers of the second kind, normalized generalized 
Bernoulli numbers B^f Cauchy numbers and normalized Cauchy numbers of the first kind Ci,„. They were introduced by James 
Gregory in 1670 in the context of area's interpolation formula (which is known nowadays as Gregory's interpolation formula) and 
were subsequently rediscovered in various contexts by many famous mathematicians, including Gregorio Fontana, Lorenzo 
Mascheroni, Pierre-Simon Lap lace, Augustin-Louis Cau chy, Jacques Binet, Ernst Schroder, Oskar Schlomilch, Charles Hermite, 
Jan C. Kluyver and Joseph Ser jl4d. vol. II, pp. 208-209], llZsl v ol. 1, p. 46, l etter written on November 23,1670 to Joh n Co llins], 


(8^ pp. 266-267, 284], 
pp. 192-194], II]i.[l8 
pm 240-251], 11^,^ 
1^ , p. 190], 0 
pp. 293-294, n° 13], li 


3. 75-78], 13S pp. 395-396], 
. 1163] . S PP- 65, 69], Slii 



11191. pp. 21-23], ln3, T. IV, pp. 205-207], ^ pp. 53-55], IlSlIl, ^ 
For more information about these important coefficients, see 1135 


p. 132, Eq. (6), m 138], |9lJ p. 258, Eq. (14 )], [93 


. 266-267, 277-280], (iMllMlp f llZll PP-106-107], 
vol. Ill, pp. 257-259 1, |167|1. llO^ p. 229], (S n“600, p. 87], 99, p. 216, n‘>75-a] 13, 
[23], Eii. Eq. (3)], |l2i|l37B pp. 128-129],Tii Chapt. 4], E3E1- 
^^These numbers, called by some authors signless generalized Bernoulli numbers \ | and signless Norlund numbers, are m uch 

less famous than Gregory's coefficients G„, but the ir study i s also very interesting, see Il35l pp. 150-151], liR p. 12], liodl , 0. 


vol. Ill, pp. 257-259], (TOj pp. 293-294, n“ 13], I811.l6l. ll92l.ll4a. ll 

^^Strictly speaking, Binet found only first four coefficients of the corresponding series for the digamma function and incor¬ 
rectly calculated the last coefficient (for K(5) he took ^ instead of ^ 1171. p. 237]), but otherwise his method and derivations 
are correct. It is also notable that Binet related coefficients K{n) to the Stirling numbers and provided two different ways for 
their computation, see Ei Einal remark]. 
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CO 2"-l „ 

^ ^ ^ ~ RkS-1 {2k + l){2k + 2) 


oo 2"+l 

= i-E E 

n=lk=2"+l 


(- 1 ) 


’^+'^n _ 1 43 

~ 12 ~ 4M 


20431 


2150797323119 


9 020112 358 835 722 225 404 403 


240240 36100888223400 236453376820564453502272 320 


( 36 ) 


was given in the first form by Niels Nielsen in 1897 11311 Eq. (6)], and in fhe second form by Ernst 
Jacobsthal in 190 6 1851. Eqs. (8)]. The same series (in various forms) was independenfly obfained by 
Addison in 1967 and by Gerst in 1969 |^]. The famous series 


T=E—= 2 

n—2 


1 1 
3 2 


2 1 
5 3 


+ ... 


(37) 


was first given by Ernst Jacobsthal in 1906 ll85l. Eqs. (9 )1 an d subsequently rediscovered by many 


1 12|] and Gerst |6 


T-n 


writers, including Giovanni Vacca llSOll . H.E. Sandham I153ll . D.E. Barrow, M.S. Klamkin, N. Miller 


Series 

OO 

7 = E 




(38) 


m — 1, n = multiple of m 

— 1, n 7^ multiple of m 

which generalizes foregoing Jacobsfhal-Vacca's series (l37t . is due to Jan C. Kluyver who discovered 
it in 1924 Q. 

In contrast, as concerns generalized Euler's constants 7m the results are much more modest. In 
1912 Hardy 1731] . by trying to generalize Jacobsthal-Vacca's series ll37ll to first generalized Euler's 
constant, obtained the following series 


1n9 °° t —It” 

7i = — £; [log2 n\ ■ (2 log2 n - [log2 2nJ) 


(39) 


n=2 


which is, however, nof a full generalization of Jacobsthal-Vacca's series since it also contains irrational 
coefficients. In 1924-1927, Kluyver tried, on several occasions ll94l. 19611 , to better Hardy's result and to 
obtain series for 7m with rational terms only, but these attempts were not successful. There also exisf 
formulae similar to Hardy's series. For instance, 

'-71+ 

= ' -72-i7^-27i7+iC'(3), 

, -73 - - 3727 - 37i 7^ + f ^(4), 

see e.g. jsZ. 120 ! ^^esides, several asymptotical representations similar to Hardy's formula are also 
known for jm- For instance, Nikolai M. Gunther and Rodion O. Kuzmin gave the following formula 


— (7 + Inn)*^ 


n—1 


m = 1 
m = 2 
m = 3 


" Ink 

71 = E-p 

k=i ^ 


In^n 


Inn 

'Yn 


1 — Inn 
12 n2 


e- 


11 — 6 Inn 
720 n4 


(40) 


^■*11 should be remarked that this series is often incorrectly attributed to Vacca, who only rediscovered it. This erro r initially 
is due to Glaisher, Hardy and Kluyver, see e.g. their works I^,l73l.l^ . It was only much later that Stefan Kramer jiosll correctly 
attributed this series to Jacobsthal i^ . 

^®Cases m = 1 and m = 2 are discussed in cited references. Formula for m = 3 was kindly communicated to the author by 
Roberto Tauraso. 
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where 0 < 0 < 1, see 1711 n° 388] o M. I. Israilov 1841] generalized expression ll40ll and showed that the 
mth Stieltjes constant may be given by a similar semi-convergent asymptotical series 


n 


7 m = E 
k=l 


In^k In’^+^n 

ln”'n 

N-l 

Blk 


{2k-l 



k m + 1 

In 

/c=l 

{2k)\ 

X 

x—n 

(2N)! 

X 


(41) 


where m = 0 , 1,2 , ..., 0 < 6 < 1, and integers n and N may be arbitrary chosen provided that N 
remains finiteUsing various series representations for the ^-function, it is also possible to obtain 
corresponding series for the Stieltjes constants. For instance, from Ser's series for the ^-function (see 
footnote 1^, it follows that 


ym — 


-1 OO -1 ?2 

n=0 k=0 


,fn yn'”+\k + l) 

kj k + 1 


m = 0,1,2. 


(42) 


An equivalent result was given by Donal F. Connon | ^ ^ 


ym — 


m 


1 00 1 n /^\ 


n—0 




(43) 


who 

asS 


used Hasse's series for the ^-function§§ Similarly, using another Ser's series expansions for 


we conclude that 


7 m — E I ^> 1+1 

n=0 




m = 0,1,2... 


(44) 


the third edition of 17111 . published in 1947, there are two errors in exercise n“ 388: the value of 71 is given as —0.073927... 
instead of —0.072815..., and the denominator of the last term has the value 176 instead of 720. These errors were corrected in 
the fourth edition of this book, published in 1951. 

^^Note that at fixed N, the greater the number n, the more accurate this formula; at n —>■ 00 it straightforwardly reduces to 
m It seems also appropriate to note here that, although Gunther, Kuzmin and Israilov obtained <40> and I41> independently, 
both these formulte ma y b e readily derived from an old semi-convergent series for the ^-function, given, for example, by 
Jorgen P. Gram in 1895 1691 p. 304, 2nd unnumbered formula] (this series for ^(s) may be, of course, much older since it is a 
simple application of the Euler-Maclaurin summation formula; note also that Gram uses a slightly different convention for the 
Bernoulli numbers). 

^®ln Eq. (3)], there is a misprint: in the denominator of the second sum 2r should be replaced by (2r ) \ [th is formula ap¬ 
pears correctly on p. 101 IsSi , but with a misprint in Eq. (3) on p. 98]. This misprint was later reproduced in Il5l]. Theorem 0.3]. 

^^Strictly speaking, Connon ItiI, p. 1] gave this formula for the generalized Stieltjes constants of which ordinary 

Stieltjes constants are simple particular cases jm = 7 m(l), see e.g. (iS Eqs. (l)-(2)]. 

^®The series representation for ^(s), which is usually attributed to Helmut Hasse, is actually due to Joseph Ser who derived 
it in 1926 in a slightly different form. The series in question is 




1 

S — 1 



{k+ir 


1 

s — 1 





The first variant was given by Ser in 1926 in jl62l. p. 1076, Eq. (7)], the second variant was given by Hasse in 1930 pp. 460- 
461]. The equivalence between two forms follows from the recurrence relation for the binomial coefficients. It is interesting that 
man y writers do not realize that Ser's formula and Hasse's formula are actually the same (the problem is also that Ser's paper 
I 162 I 1 contains errors, e.g. in Eq. (2), p. 1075, the last term in the second line should be ( —l)"(n + 1 )^’’ i nstead of (—l)"n^®). An 
equivalent series r epre sentation was also independently discovered by Jonathan Sondow in 1994 Il6,5ll . 

^^Ser's formula 1162], p. 1076, Eq. (4)], corrected (see footnote|30) and written in our notations, reads 


as) = 


- |G,j+i| X^(-l)* 

n=0 A:=0 


(k + l)- 
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where G„ are Gregory's coefficients, see (l3TTl - ll33tl^ This series was also independently discovered 
by Marc-Antoine Coppo in 1997 11451, p. 355, Eq. (5)] by the method of finite differences. Several m ore 
complicated series representations for jm with irrational coefficients may be found in [iziEiiS 


168,128, 


186, 


38, 


390 


111.2. Derivation of the series expansion 

Consider now series l fT9ll . A formal rearrangement of this expression may produce a series with 
rational terms only for 7^- In view of the fact that 


m+1) = E 


n—k 


I Si (n. A:) I 


n ■ n\ 


= E 


I Si (n. A:) I 


n—1 


see e.g. |92, pp. 166,194-195]and that 


n ■ n\ 


fc = 1,2,3,... 


^ ^ ^ l-{lk)\ 2-{lk)\ ' 

the formal interchanging of the order of summation in lfT9t leads to 


k = 1,2,3,... 


(45) 


(-l)'”m! g 1 


(-l)fc|Si(2fc + 2,m + l)| ■ |Si(n,2fc + l)| 


n 


n-n\ 

n=l k=0 


(27t)2^+i 


(-l)™m! ^ (-l)^|Si(2A + 2,m + l)| 


TC 


k=0 


= (-ir'^iE' 


(27T)2*^+1 

Si (2A!, T 1) I ■ 


E 

n—1 


|Si(n,2A:+l)| 


n ■ n\ 


aik+i) 


{2k)\ 


Remarking that such operations are not permitted when series are not absolutely convergent (that is 
why we wrote x instead of =), we imderstand why the resulting series diverges. However, since 


^^Note that for m = 0, the latter series reduces to Fontana-Mascheroni's series im. ,_, 

Apart from formula lUl . in Israilov also obtained several other series representations for 7 ^. Stankus llbsll showed 
that first two Stieltjes constants may be represented by the series involving the divisor function. Works of several au thor s 
showing that Stieltjes constants are related to series containing nontrivial zeros of the ^-function are summarized in lisdi . 
Coffey gave various series representations for 7 ^ in However, we also noted that 1^.1^ both contain numerous 

rediscoveries, as well as inaccuracies in attribution of formulee (see also liiil l. For instance, "Addison-type series" Eq. (1.3)] 
are actually due to Nielsen and Jacobsthal, see our Eq. <36t above. Numbers Pn+i are si mply signless Gregory's coeffi cient s 
IGnl and their asymptotics, Eq. (4.8)/(2.82) (ii, p. 473/31], is known at least since 1879 Il6l]. Eqs. (25)-(25a)] (see also Il70ll . 
Il7l]. pp. 106-107] ). Th eir "full" asymptotics, Eq. (4.^/(2.84) p. 473/31], is also known and was given, for example, by 
Van Veen in 1951 II81I1 , lisdl . Representation (2.17) Ijg, p. 455/13Hs due to Hermite, see Eq. (13) ll^. p. 5 41]. Formula (1.17) 
13^, p . 2052] was already obtained at least twice: by Binet in 1839 Iia p. 257] and by Norlund in 1923 Il35l. p. 244], see details 
in list Final remarkl. Formula (1.18) Ij^. p. 2052] straightforwardly follows from Ser's results I162I] dating back to 1926, etc. 

^^See also Il63l,ll54ll , where this important result was rediscovered much later. By the way, this formula may be generalized 
to the Hurwitz Gfnnction 

(^{k + l,v) = r 1^1 (”'m i: = 1 , 2 , 3 ,... , Rez?>0, 

see Gi] . At n —t 00 , the general term of this series is O 



16 






























this series is alternating, for any prescribed m and N, one can always find such 9 G (0,1), generally 
depending on m and N, thaf 


N 


Tm — 2^m,0 


{-irml-i: 

k=l 


I Si {2k, m + 1) I ■ B2k 


+ 


e- 


|Si(2N + 2,m + l)| ■B2N+2 

{2N + 2)\ 


1 

12 


1 , 1 


1 , 1 


11 


137 


+ ■ 


121 


7129 


57844301 


= 


12 720 15120 11200 332 640 908107200 

469 6515 


+o-l + A_ 

60 336 21600 


131 672 123 63 427 

+ 


-0 


133056 825552 000 89100 

1 17 , 967 4523 , 33 735311 9301169 

120 ~ 1008 28800 ~ 49896 101088000 ~ 5 702400 


m = 0 


m = 1 


m = 2 


m = 3 


(46) 


holds sfrictlyo Moreover, taking into account dl^ , above series may be always written in a form 
wifhout Stirling numbers. For instance, for Euler's constant and for first three Stieltjes constants, it 
becomes 


N 


i = +lh+E-¥+»- 


k=l 


1 ^ B2fc ■ H2k-1 , , 
T'l = -sL - Z - 


B 2 N +2 [ 

N + 1 j 

B2N+2 ■ H2N+1 


■ k=l 


2N + 2 


1^ ^2k' {^2k-l ^2J:-l} , a ®2N+2 ' {WIn+I 

■5 - 1 - + ^- 


72 = 


(47) 


■ k=l 


2N + 2 


1 N B,, . - 3H2,.,Hfl, + 2Hfl,) 

= “2 ^- k - 

^ k=l ^ 


■9- 


B 2 N +2 ■ - 3H2n+iH§^, + 2H^"4J 


A3) 


2N + 2 

where 0 < 0 < 1 and A1 < 00 (these parameters are different in all equations, and in each equation 
9, in general, depends on N). By the way, one may notice that first of these formulae coincide with 
Stirling series ll29ll . while other formulae are, to our knowledge, new and seem to be never released 
before. Derived formal series are alternating and are sometimes referred to as semi-convergent series or 


^®There is another way to obtain (46). Consider first till at n = 1, and then use 1131 to show that 
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Figure 3: Partial sums of series f47l for 7 and 71 (on the left and on the right respectively) for N = 1,2,..., 10. 
Blue lines indicate the true value of 7 and 71 , while dashed black lines with the red points display corresponding 
partial sums given by i47l . 


divergent enveloping series 


E 


It is also easy to see that they diverge very rapidly 


Si {2k, m + 1) I ■ B2k 


,_i (2fc-l)!.ln"(2fc-l) 
ml ■ ( 27 r)^^ 


2^/W{-lf-^ IrTk 
ml VF 


-) 

Tie J 


Ik 


fc —t 00, m = 0,1,2,..., so rapidly that even the corresponding power series Y, 


Si(2*:,m+1) 
{2k)\ 


B2kX 


.2k 


diverges everywhere!^ Behaviour of this series for first two Stieltjes constants is shown in Fig.|3] 


111.3. Some series transformations applied to the derived divergent series 

In order to convert ll46t into a convergent series, one may try to apply various series transforma¬ 
tions and regularization procedures. However, since (l46l is strongly divergent, the use of standard 
summation methods may not result in a convergent series. For example, applying Euler's transfor- 
matiorj^ we obtain another series with rational coefficients only 


7m — + { — ly 


1 

k=l ^ n=l 


N 

"ml^ 


|Si(2n, ?n -h 1) I ■ B2n 


k-1 
n — 1 


+ R,n{N) = 


^^These series were an object of study of almost all grea t mathematicians; t he re ader interested in a deeper study of these 
series may wish to consult the followi ng lit erature: I22l.l5ll. l^. 1^. Chapt. XI], jl4ll. Chapt. 4, §1], exercises n° 374-n° 388 
in El PP- 46-48], ll^.[i3i[i^. H.l^lm.y. [i43.[]5 ]^.lM 


^^Coefficients jSi {2k, m + 1) | and Bernoulli numbers both grow ve ry qu ickly: as t 00 we haye [Si {2k, m + 1) | ~ {2k — 
1)! In’"]!*: - h/m!, see E), and 624 ~ 2(-l)*’-i(27r)-2^)!, see e.g. p. 5], EE P- 261]. 

38See e.g. S pp. 244-246], pp. 144 & 170-171], ES pp. 269-278 & 305-306]. 
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Figure 4: Euler's transformation of series f47t for 7 and 71 (on the left and on the right respectively) for N = 
1,2,..., 14. Blue lines indicate the true value of 7 and 71 , while dashed black lines with the red points display 
corresponding partial sums given by (l48t . 


f 1 , 1 , 3 , 89 , 37 , 299 

2 ^ 1^ 10080 8960 147840 ' 


1 49 187 5431 91151 


24 

2880 

24192 1612 800 53222400 

1 

31 

4093 

50 789 

602325403 

240 

13440 

2419200 

106444800 

581188608000 

1 

1 

1609 

120749 

694 773847 

480 

40320 

3225600 

159667200 

498161664000 


m = 0 

m = 1 

(48) 

m = 2 

ni = 3 


which are all divergent (their remainder R(N) —t 00 as N —t 00). At the same time, these series behave 
much better than (l46t . Thus, the series for 7 starts to clearly diverge only from N ^ 10, and that for 
7i only from N ^ 8, see Fig. S) The minimum error for the first series corresponds to N = 7 and 
equals 3 x 10“^, that for 71 also corresponds to M = 7 and equals 9 x 10“^. Attempts to regularize 
series ll46t - ll47t with the help of Cesaro summation are also fruitless since its general term grows more 
rapidly than fc at A: —t 00. Similarly, BoreTs summation does not provide a convergent result. 


777.4. An estimate for generalized Euler's constants 

Finally, we note that derived formal series dlbll provides an estimation for generalized Euler's 
constants. Since this series is enveloping, its true value always lies between two neighbouring partial 
sums. If, for example, we retain only the first non-vanishing term, the mth Stieltjes constant will 
be stretched between the first and second non-vanishing partial sums. Thus, accounting for known 
properties of the Stirling numbers of the first kind 

|Si(m -F 1, m -F 1) I = 1 

|Si(m-F2,m-F 1)| = 2{m + l){m + 2) 

|Si(m-F3, m-F 1)| = ^(m-F l)(m -F 2)(m -F 3)(3ot-F 8) 
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|Si(m + 4, m + 1) I = + l)(m + 2)(m + 3)^(m + 4)^ 

which are valid for m = 1,2,..., we have 


^m+l 


m • 




(3m+ 8) ■ \ Bm+3\ 

24 


B 


m+l 


I ®m+i I (3m + 8) ■ I Bffi^31 


m + l 
Bm+2 


24 




m +1 ' 

Bm+1 


m +1 ' 




(m + 3) (m + 4) ■ 16^+41 


B 


m+2 


2 ' ' 48 

|Bm+2| (m + 3){m + A) ■ \Bm+i\ 


+ 'Yin + 


2 

I Bm+21 


m = 1,5,9,... 


m = 3,7,11,... 


m = 2,6, 10 ,... 


m = 4,8,12,... 


( 49 ) 


2 48 ' 2 

Case m = 4,8,12,... may be also extended to m = 0, if we recall that in this case it gives bounds 
for 'Y — 2 , see This case yields < 7 < ^, which is undoubtedly true. Note also that above 
bounds are always rational, which may be of interest m certain circumstances. 

Estimation (|4^ is relatively tight for moderate values of m, and becomes less and less accurate 
as m increases. However, even for large m, it remains more accurate than the well-known Berndt's 
estimation 

2(m- 1)! 


7ni ^ 


T^m 

4 (m — 1)! 

TT"* 


m = 1,3,5,... 


m = 2,4,6,... 


(50) 


see Bid pp. 152-153], more accurate than Lavrik's estimation |7m| ^ see 111131 Lemma 4], 

more accurate than Israilov's estimation 

m! C{k) 


ITj) 


< 


{2kY 


for k = 1,2,3, where C(l) = j, C(2) = C(3) = see 

Williams' estimation 

2 (2m)! 


\'Yni\ ^ 


m 


m+l 


m 


{2nY 
4 (2m)! 
{2nY 


(51) 

, and more accurate than Nan-You- 
m = 1,3,5,... 


(52) 


m+l 


m = 2,4 ,6,... 


see Bl28l. pp. 148-149]. Besides, our estimation also contains a si gn, while above estimations are sign¬ 
less. At the same time. (l4^ is worse than Matsuoka's estimation lll2lLll22 l. |7m| < 10“^ln™m,m ^ 5, 
which, as far as we know, is currently the best known estimation in terms of elementary functions for 
the Stieltjes constants 1^ Note, however, that estimation's bounds ll49t may be bettered if we transform 
parent series (l46l into a less divergent series, provided the new series remains enveloping. 


^^Numerical simulations suggest that Matsuoka's estimation I 121 I.I 122 I 1 may be considerably improved, see e.g. Re¬ 
cently, Knessl and Coffey reported that they succeeded to significantly better Matsuoka's estimation and even to predict the 
sign of 7m. The authors published their findings in j^, and also reprinted them in Theorem 1]. We, however, were not 
able to verify these results, because several important details related t o v( n) from pp. 1 79-180 I 9 I 1 were omitted. Estimation 
of the similar nature was later proposed by Adell j^I, but Saad Eddin jl5ol. Tab. 2], jlSlIl reported that AdelTs estimation may 
provide less a ccura te results than Matsuoka's estimation. Saad Eddin also provides an interesting estimation for the Stieltjes 
constants, see IlSOl, Tab. 2] and mentions some further works related to the estimations of the derivatives of certain L-functions. 
Yet, very recently we foimd another work devoted to the estimation of Stieltjes constants 1^ : the latter resorts to the Lambert 
W-function. 
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Appendix A. Two simple integral formulae for the Stirling numbers of the first kind 


Appendix A. 1. The first formula 

Consider equation ([T^) defining signless Stirling numbers of the first kind. Dividing both sides by 
where k = 1,2,3,, and integrating along a simple closed curve L encircling the origin in the 
counterclockwise direction, we have 



L 


dz 


^ |Si(n,/)| ■ ^ ^ dz = 2m\Si{n,k)\ 

/=i 


Ini Si t 


(A.l) 


in virtue of Cauchy's theorem. In practice, it is common to take as L the unit circle, and hence 


Si{n,k) 


Ini 


ft + TT 

T z^+i 2n J r(e'‘P) ^ 

DL—n 


|z|=l 


(A.2) 


where a e R. 

Alternatively, in virtue of the Cauchy residue theorem, we also have 

|Si(n,fc)| = res , n,k = 1,2,3,... (A.3) 

' ' z=o z'^+^r(z) 


Appendix A.2. The second formida 

Consider generating equation (|8h) for the unsigned Stirling numbers of the first kind. Proceeding 
as above and then making a change of variable z = re"?', we obtain 


|Si(n,A:)| = 


(-1)^ 

2Tzi 


n\ 

M 


In^il-z) 


yir + l 


dz = 


, tt + TT , , . ^ 

(-1)^ ^ f \n\l-re^V) 


i(pn 


2n k\ 


yfl 


dcp (A.4) 


|z|=r 


where, due to the radius of convergence of 1(8^), 0 < r < 1, and a € R. Obviously, the line integral 
in the middle may be taken not only along the indicated circle, but along any simple closed curve 
encircling the origin in the counterclockwise direction and lying inside the unit circle |z| = 1. By the 
way, it is interesting that the same integral taken between (0,1) reduces to a finite linear combination 
of Stirling numbers of the first kind and ^-functions, see Sect. 2.2, Eq. (48)]. It seems also ap¬ 
propriate to note here that several slightly different integral formulae of the same kind as llA.4b were 
given by Dingle lISOl. pp. 92, 199]. The author, howe ver, d id not specify the integration contour. This 
inaccuracy has been partially corrected by Temme 11174 p. 237], who indicated that the integration 
contour should be a "small circle"!^ 

Similarly to the previous case, the signless Stirling numbers of the first kind may be also given by 
the following residue 




n,k = 1,2,3,.. 


(A.5) 


^®We, however, note this condition is not really necessary; it is sufficient that the integration path be a simple closed curve 
encircling the origin in the right direction and lying inside the unit disc. This remark may be important for the numerical 
evaluation of these integrals, which have been reported as not well-suited for these purposes. 
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Figure A.5: Relative error of the estimation for the Stirling numbers of the first kind given by <A.6t as a function 
of k for n = 80, logarithmic scale. 


Finally, expression (IA.4b may be also used to readily get an interesting botmd for Stirling num¬ 
bers. In fact, one may notice that if r ^ 1 — e~^ w 0.63, then the principal branch | ln^(l — z)| ^ 1 
independently of k. Therefore, for k = 1,2,...,n 


Si(M,fc)| < 



\z\=r 


ln^(l — z) 
2« + l 


dz\ ^ 


nl 

(l-e-i)”k! 


(A.6) 


since |dz| = rdcp. If n is prescribed, this estimation is relatively rough for small and large (close to 
n) factors k; in contrast, for values of k which are slightly greater than n/2, this estimation is quite 
accurate, see Fig. IA.5l 
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